Resonant excitation of graphene K-phonon and intra-Landau level excitons in 

magneto-optical spectroscopy. 
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Precise infra-red magneto-transmission experiments have been performed in magnetic fields up 
to 32 T on a series of multi-layer epitaxial graphene samples. We observe changes in the spectral 
features and broadening of the main cyclotron transition when the incoming photon energy is in 
resonance with the lowest Landau level separation and the energy of a K point optical phonon. We 
have developed a theory that explains and quantitatively reproduces the frequency and magnetic 
field dependence of the phenomenon as the absorption of a photon together with the simultaneous 
creation of an intervalley, intra-Landau level exciton and a if-phonon. 



Phonons, electrons, holes and plasmons are fundamen- 
tal elementary excitations in solids, and their interactions 
among themselves and with light strongly influence the 
electronic and optical properties of a system. Graphene, 
a single atomic layer of carbon, exhibits many fascinat- 
ing properties of basic and practical interest. Its low- 
energy charge carriers behave like two-dimensional (2D) 
massless Dirac fermions. The unique electronic struc- 
ture of graphene has led to the prediction and obser- 
vation of a number of novel phenomena involving var- 
ious multi-elementary excitation interactions including 
electron-photon, electron-phonon, electron-electron, and 
electron-plasmon couplings. Recent discoveries include 
the observation of plasmarons in ARPES |l|, controlling 
of quantum pathways in inelastic light scattering Q , and 
strong excitonic effects in optical absorption , among 
others. 

In this paper, we report the discovery of another in- 
triguing phenomenon seen in the magneto-optical re- 
sponse of graphene when the incoming photon energy is 
in resonance with the lowest Landau level separation and 
the energy of a Brillouin-zone-boundary phonon at the 
K point. The frequency and magnetic field dependence 
of the spectral features are understood and explained in 
terms of resonant transitions which involve the absorp- 
tion of a photon together with the simultaneous creation 
of an intervalley, intra-Landau level exciton and a K- 
phonon. The phenomenon provides a novel manifestation 
of mult i- elementary excitation interactions in condensed 
matter. 

The phonons that interact the most strongly with 
electrons in graphene are those at the Brillouin-zone- 
center (L) and at the Brillouin-zone-boundary (A) [6(. 
The T-phonon has been the subject of theoretical study 
Q and its signatures have been observed in several 



experiments 0, Q. The K- phonon also interacts strongly 
with electrons, giving rise to characteristic features in 
Raman spectra [8(. Compared to L-phonons, the effect 
of A-phonons on optical experiments is more subtle and 
indirect, due to the large difference in momenta between 
photons and the A-phonons. The large wavevector of 
the A-phonons may also couple the two inequivalent val- 
leys in the graphene band structure, introducing an ad- 
ditional degree of freedom into the system. The present 
study shows that this valley degree of freedom in the 
coupling has unique and important consequences in the 
magneto-optical response of graphene. 

In our experiment, precise infra-red transmission 
measurements were performed on multi-layer epitaxial 
graphene samples, at 1.8 K, under magnetic fields up to 
32 T. It is known that under a perpendicular magnetic 
field B, the electronic structure of graphene is quan- 
tized to discrete Landau levels (LL) with the energy 
E n = sgn(n)vFy2ehB\n\ where n are integers including 
0. Multiple optical transitions are allowed and have been 
observed [§] between LL \n\ to \n\ ± 1 if these levels have 
appropriate occupation factors. We focus our attention 
on the magnetic field evolution of the transmission spec- 
trum near the main optical absorption peak, with energy 
Eqi = vp\/2ehB involving the n = LL (i.e., transitions 
from n = -1 to n = or from n = to n — 1.) It is 
discovered that a strong change in the optical spectrum 
occurs when Eqi reaches an energy near the energy of a 
A"-phonon. This is interpreted as a signature of a new 
multi-elementary excitation phenomenon involving inter- 
action between electrons and the A-phonon, and a theory 
is developed to explain quantitatively the phenomenon. 

A series of four multi-layer epitaxial graphene sam- 
ples were grown [To| on the C-terminated surface of SiC. 
The thickness of the SiC substrate has been reduced 
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FIG. 1: (a) Relative transmission spectra (T(B)/T(0)) of sample SI, for different magnetic field values, (b) Same data as in 
(a) presented in a false color plot of -ln[T(B)/T(0)]. (c) Bottom panel: variation of cyclotron transition energy Eoi/B 1 ^ 2 as a 
function of B 1 / 2 (open blue circles). Top panel: variation of the fitted linewidth Fio(-B) of the transition Eoi(B) as a function 
of B 1 / 2 (open red circles). This variation is de-convoluted is two parts, one named 701(B) (green dots) and the remaining part 
(stars) which is assigned to the if-phonon contribution. 



to 60/zto in order to minimize the very strong double- 
phonon absorption of SiC in the energy range of inter- 
est. Whereas the transmission spectra reveal some contri- 
butions of graphite Bernal-stacked inclusions, the major 
part displays the electronic band structure of an isolated 
graphene monolayer resulting from the characteristic ro- 
tational stacking of the sheets [Tljj . Using standard tech- 
niques [lH , we measure the relative transmission spectra 
TR(£?, uj)(see Supplementary Information SI-I). 

The TR spectra of the sample SI are displayed in 
Fig. [lj, and b for different values of B. They show, as for 
all samples, an apparent splitting of the transmission dip 
into two and a pronounced increase of the broadening of 
the dip, as well as a change in the f?-field dependence of 
the overall position Eq\ for fields larger than 17 T. 

In general, the detailed analysis of such spectra re- 
quires the use of a multi-layer dielectric model including 
all layer dielectric properties of the sample. In particular, 
for each graphene sheet, one has to introduce the corre- 
sponding components of the optical conductivity tensor 
a xx (uj) and a xy (uj), defined by J a = J2p a ap£i3 (J a an d 
£p denote current density and electric field vectors re- 
spectively). Here, the x and y-axis lie in the plane of the 
sample. For instance a xx (uj), for transitions involving the 
n=0 LL, is written as: 



<7 XX (UJ,B) = I 



e 3 B^ 



MUU{B) - fs(B)) 
huj ^—^ hw — E ryS (B) + iT rs {B) 



(1) 



where r,s scan the values and ±1, < f r < 1 is the 
occupation factor of the LL r, M TjS the optical matrix 
element, E Tj8 = E r — E s = E01 and T rs (B) — r i(.B) 
measures the broadening of the transition . However this 
approach requires the knowledge of the number of effec- 



tive active layers as well as their carrier densities which, 
in turn, implies some approximations (see SI-II for de- 
tails). We will use this approach in a refined analysis, 
but in our initial analysis, fit directly the Eqi transition 
with a single Lorentzian line. 

This fit provides, for each value of B, two indepen- 
dent parameters Eqi(B) and Foi(-B) which are plotted 
in Fig. [TJ: as a function of B 1 !" 1 . The variation of Foi(-B) 
shows clearly an extra bump beyond 17 T (clearly ap- 
parent in Fig. [lb ) whereas that of E m (B)/ B 1 / 2 ex- 
hibits a downwards kink at the same field. The evolu- 
tion of Tqi(B) is decomposed in two parts: one named 
701 (meV) = 2.3 + 0.6-/B + 0.245 (full dots) and an ex- 
tra contribution to r i represented by stars in the top 
panel of Fig. [TJ:. The same decomposition is obtained 
for all samples with slightly different coefficients which 
are discussed in SI-II. The \f~B~ dependence of 701 is at- 
tributed to scattering with short range impurities [l4j : 
the linear variation with B is not predicted by scattering 
mechanisms considered in previous studies 14| and will 
be discussed later. The shape of the extra contribution 
to Toi has a threshold at around 17 T (~ 151 meV) with 
a maximum amplitude of about 2 meV for all samples, 
within the experimental uncertainties, before decreasing 
at higher fields. In that energy range the only likely 
excitations in graphene which could play a role are the 
zone-boundary if-phonons 15, Hiij . Note that, at this 
threshold energy, the evolution Eqi/B 1 / 2 (bottom panel 
of Fig. [TJ:) is the signature of interactions because, in 
the absence of such interactions, it should be an horizon- 
tal line. This is reminiscent of the response function of 
the Frohlich interaction observed in polar semiconductors 
12 . 13| where a zone-center longitudinal optical phonon 
is emitted when the energy of the cyclotron transition 
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FIG. 2: Cyclotron resonant excitations \r)i(j, k))(a) and 
\r)2{j, k)) (b) involving the n = Landau levels, at either 
the K or if' valley, with the excitation energy E ex = Ecu. (c) 
Excited state |f<f(,j, fe)) corresponding to intervalley transfer of 
an electron within the n — LL from .K^ to 7?' and creation of 
a phonon with wavevector K + q, which has excitation energy 
E ex = fihJg+g ■ (d) At resonance, i.e. Eoi = fkug + ^, \rji) and 

\r]2) may be transformed to |£) via an emission of a K + q 
phonon as indicated by the dashed brown lines. The mixing 
between the three excited states, mediated by the electron- 
phonon interaction, is the strongest when the resonance condi- 
tion is satisfied (-Eoi = ^k-h?) • Electron- phonon scatterings 
involving LL |ra| > 2 do not result in similar resonances in the 
experimental energy range (see SI-IV.) In general, there are 
scattering processes involving phonons of different q vectors 
near K, but we find the coupling strength significant only in 
a small region with q « where the phonon energy is essen- 
tially constant and equal to Tiwk- The theoretical description 
of this continuum of near K phonons may be incorporated 
effectively by considering a single state |£) with a complex 
self energy and using an appropriately renormalized coupling 
constant between |?7i), |jj2) and \Q (see derivation in SI-IV) . 



exceeds that of the phonon. However, graphene is not 
a polar material and, furthermore, we are not expecting 
a direct interaction of the A"-phonon with the infra-red 
light. One therefore has to explain why and how the 
if-phonons enter the interaction. 

The mechanism for the interaction of the cyclotron 
transition with the if -phonon is not a priori clear for 
the reason of conservation of wavevector, i.e., a state cor- 
responding to an optical cyclotron transition which has 
wavevector cannot decay into a state with a phonon of 
wavevector K without simultaneously emitting another 
entity of wavevector —K. Given that this interaction 
must necessarily involve the electron-phonon interaction, 
and that the low-energy electronic states at the two val- 
leys of graphene are separated in reciprocal space by 
a wavevector if, we propose that an intervalley intra- 



LL electron scattering process (creation of a zero-energy 
intervalley electron-hole pair) is responsible for the ob- 
served phenomenon (Fig. [5J) 

To explain the observed spectral features, we shall con- 
sider the coupling of the photo-excited (cyclotron tran- 
sition) states (denoted as \rji(j,k)) and \r]2(j,k)) for the 
— > 1 and — 1 — > LL transitions respectively) with 
an excited state of the system containing an intervalley, 
intra-Landau level electron-hole pair excitation and a if- 
phonon of wavevector K + q (denoted as \£$(j, k))). See 
Fig. 2. Here j is the valley index and k is the quantum 
number which describes the degenerate states within a 
LL. As illustrated in Fig. 2, at cyclotron transition energy 
-Eoi near huig + ~, \r)i(j, k)} and 1 972 ( J , k)} are connected to 

|£<f(i> k)) via an emission of if +q phonons. The electron- 
phonon coupling in the system can then strongly mix 
these states at resonance and lead to significant changes 
in the absorption spectrum. 

In order to reproduce quantitatively the transmission 
spectra, one has now to use the multi-layer dielectric 
model (see SI-II for details). The optical conductivity 
is calculated using the Green's function formalism intro- 
duced by Toyozawa (l8| . The conductivity components 
are 



Rea tj (huj) oc Im(0|M i G(/kj)M : ,|0) 



(2) 



where |0) is the ground state of the system, and Mj, 
Mj are row and column vectors containing optical ma- 
trix elements. G is the retarded Green's function written 
here as G(hw) = [tkJ — H — E(w) + is}^ 1 with s — > 0~ 
with H being a 3x3 Hamiltonian describing the inter- 
action between the three states \r] 2 ) and The 
poles of G are broadened by a self energy reflecting 
additional lifetime effects due to the if-phonon disper- 
sion and the excitations from electron-acoustic phonons 
scattering, electron-electron interactions and other envi- 
ronmental effects not captured by H (see SI-V). In the 
theoretical calculations, we have incorporated this addi- 
tional physics near resonance with a broadening func- 
tion of the form Tqi(lj,B) = 701(E) + Im(S(cj)) with 
701(E) given above and Im(E) taken to be Im(S(w)) — 
TrA 2 f(R/fuv K )9(huj - huj K )e- R{ul/uiK - 1 '> . The derivation 
of the form of Im (S(w)) within our theoretical model is 
given in SI-IV and SI-V. A 2 ( meV 2 ) = 1.94 x E(T) is the 
electron-phonon parameter determined from Ref. [l6j |. 
/ is a factor, dependent on the relative magnitudes of 
the valley (Ay) and spin (As) spittings of the n = LL, 
that is determined by the occupations of the spin and val- 
ley sublevels of the n = LL (see SI-IV). The threshold 
structure of Im(E(w)) represents the onset of if -phonon 
emission at energies above hwx- The only parameters in 
the theory are the phonon energy Hlok and the factor R. 
The absorption spectra can be obtained from Re a(u>) by 
a simple proportionality relation (see SI-V and j§|). 

In Fig. [3] we compare the experimental data, for sam- 
ple SI, with the calculated results from the model with 
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FIG. 3: Comparison, for different values of the magnetic field, 
of experimental data for sample SI (open black dots) with the 
theoretical curves (red lines) using the proposed K-phonon 
interaction model. As a guide for the eye, upward arrows show 
the evolution of the main cyclotron resonance line whereas 
downward arrows point the structure corresponding to the 
emission of A-phonons. 



Ay > Ag. We have taken Hu>k = 151 meV and R = 3. 
The Fermi velocity vfoi = 1.012 x 10 6 ms _1 for the tran- 
sitions between the n = LL and n = ±1 LL has been 
determined experimentally from the positions of the cy- 
clotron transition lines (see SFII and SITU.) The same 
parameters have been used for all the graphene samples 
(see SITU) . The good agreement between theory and ex- 
periment, for all samples, lends strong support to the 
physical mechanism proposed for the phenomenon. It 
is , however only obtained when assuming Ay > As 
(Fig.©. 

An interpretation of the linear term in 701 (B) is that 
it may be a broadening resulting from the breaking of 
the valley degeneracy (Fig. consistent with the pre- 
diction that Ay oc B, according to several theoretical 
models (2p| - [22l |. In this picture, the values of Eqi at the 
two valleys differ by Ay 



23j, resulting in a broadening 



of the main line by a similar amount. An inspection of 
the linear term in 701 (B) shows that it is greater than 
the value of Ag due to Zeeman splitting. This is consis- 
tent with our conclusion that Ay > A5 in our samples 
(Fig. [4|. However, further investigations are needed to 
test if the valley splitting is indeed the main contributor 
to the linear term in 701. 

To conclude, we have observed splitting of spectral fea- 
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FIG. 4: (a), (b) The spin and valley sublevels of the n = LL, 
for the case when Av > As and Ay < As respectively. In 
the former case, intervalley electronic transitions are possible 
for both spin orientations. In the latter case, only one spin ori- 
entation supports intervalley electronic transitions, owing to 
the Fermi factor. The effective coupling strength is therefore 
different in the two cases giving rise to different predictions 
for the transmission spectra, (c) A comparison between ex- 
perimental spectra and theoretical simulations for Av > As 
and Av < As for values of magnetic field B = 18. 5T and 
B — 23. 0T respectively. For all our samples the splitting of 
the valley degeneracy is larger than the spin splitting. This is 
consistent with recent tunnelling spectroscopy experiments on 
graphene on SiC Theoretical discussions on the relative 
magnitudes of these splittings can be found in Refs. [20l425| ]. 



tures and broadening of the main absorption line above 
17 T in our infra-red transmission experiments on epi- 
taxial graphene that we attribute to a resonance be- 
tween cyclotron transitions and K-phonon emission cou- 
pled with a zero-energy intervalley electron-hole excita- 
tion. We have developed a theoretical model that repro- 
duces the experimental results with quantitative agree- 
ment, thereby shedding light on the nature of the optical 
transitions and electron-phonon interaction in graphene 
under strong magnetic fields. Our results indicate that 
the valley splitting of the LL is larger than the spin split- 
ting in our samples. To our knowledge, this is the first 
instance where information on the sublevel structure of 
graphene has been obtained using infra-red transmission 
measurements, c omp lementing transport- 26], [27| and 
tunnelling-based [19( experiments. 
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I. EXPERIMENTAL DETAILS 



Far infra-red magneto-optical experiments have been performed at magnetic field strengths 
up to 29-32 T and at a fixed temperature of 1.8 K. The light (provided and analyzed by a 
Fourier transform spectrometer) was delivered to the sample by means of light-pipe optics. 
All experiments were performed with nonpolarized light, in the Faraday geometry with the 
wave vector of the incoming light parallel to the magnetic field direction and perpendicular 
to the plane of the samples. A Si bolometer was placed directly beneath the sample to 
detect the transmitted radiation. The response of this bolometer is strongly dependent on 
the magnetic field. Therefore, in order to measure the absolute transmission TA (B,u>), we 
used a sample-rotating holder and measure for each value of B a reference spectrum through 
a hole. These spectra are normalized in turn with respect to TA(0, oS) to obtain a relative 
transmission spectrum TK(B,u) which only displays the magnetic field dependent features. 
Those spectra are presented in Fig. [TJ 



II. GLOBAL ANALYSIS OF THE SPECTRA 

In Fig.[U the relative transmission spectra are shown for sample S2 (a), sample S3 (b) and 
sample S4 (c). In general, for filling factors v < 2, the oscillator strength of a magneto-optical 
transition in graphene is expected to be proportional to \/B if the broadening parameter r i 
is a constant. If r i ~ y/B the minimum of the transmission should be independent on B. 
However and as for the sample SI (Fig. []Ja) of the main text), all the observed transmission 
data display an extra broadening of the Eqi transition for magnetic fields values larger than 
17 T. This leads us to conclude that such an effect is an intrinsic property of the graphene 
sheets. 
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FIG. 1. Relative transmission spectra, as a function of the energy, for different values of the 
magnetic field B : (a) sample S2, (b) sample S3 and (c) Sample S4. 
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FIG. 2. Bottom panels: variation of Ebi/V-B as a function of \/B (open blue circles). Top panels: 
variation of the fitted linewidth T\q(B) of the transition Eq\(B) as a function of \f~B (open red cir- 
cles). This variation is decomposed is two parts, one named 701 (B) (green dots) and the remaining 
part ( wine stars) which will be assigned to the .fT-phonon contribution, (a) sample S2, (b) sample 
S3 and (c) sample S4. 



One can, as said in the main paper, analyze globally the E'en transition by a fit with a 
single Lorentzian line which provides the evolution of the two parameters Eqi (B) and Toi (B) 
as a function of B 1 ^ 2 as displayed in Fig.[2]for samples S2, S3 and S4. However, as this will be 
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explained in the following, the model for the electron- .fT-phonon interaction is dependent on 
the filling factor v of the graphene layer and therefore on its carrier density Ns [y = Ns§o/B, 
$ being the flux quantum). We have therefore used a multi-layer dielectric model to analyze 
the data. This model assumes that each graphene sheet is uniformly spread over the sample. 
This is a strong assumption, difficult to justify a priori but which is necessary here to get 
some relative quantitative information. If one knows, for a given sample, the number N e ff 
of effective layers and their carrier densities, one can fit quantitatively the data using Eq.l of 
the main text when Vq, the Fermi velocity entering the optical matrix element (M rs oc Vq), 
is known. Very often the value vq is taken identical to the Fermi velocity Vp deduced from 
the magnetic field dependence of the optical transitions. However vp is the re- normalized 
velocity including electron-electron interactions and it has been shown l| that Vq should be 
the band velocity ignoring these interactions as given, for instance, by LDA calculations. In 
the present work, we have taken for all samples Vq = 0.85 x 10 6 m/sec. We next evaluate 
the number N e ff for each sample. In the range of magnetic fields 13 to 17 T, the relative 
transmission spectra (Fig. [1]) reaches values above 1 which depends on the number N e ff. 
we have therefore a guide to estimate this quantity. This is however only a guide because 
this value is also sensitive to the effective coverage of each layer over the sample and, of 
course, to experimental uncertainties. We estimate N e ff = 3 for sample S2, N e ff = 5 for 
sample S3 and N e ff = 6 for samples SI and S4. The carrier density N s for each layer is 
determined in the following way: one knows that, for 2 < v < 6, upon increasing B, the 
intensity of the Eqi absorption starts to increase, at the expense of the intensity of the E 12 
transition (E12 = E2 — E{). The intensity does not change with B for v < 2. Therefore, the 
disappearance of the optical transition E\i corresponds to v — 2. Following step by step, 
as a function of B, the transmission spectra one can evaluate the carrier density iV sm for 



each layer m. This is a loop process which converges reasonably (within 20 per cent) but 
has to be done independently for each sample. This gives reliable values for the uppermost 
layers which are the less doped and will play a major role when the electron- i^-phonon 
interaction will be switched on. The value of N a \ for the layer close to the SiC substrate 
can be set arbitrary to 5 to 6 10 12 cm -2 as given by transport data on samples grown under 
similar conditions: this layer indeed does not contribute to the transition Eqi in the present 
experiment. Fitting the data with such a multi-layer dielectric model and using Eq.l of the 
main text provides also the two parameters E m (B) and Tqi(B): it turns out that within the 
experimental errors (error bars in Fig. [2]) they are identical to the results obtained previously 
by fitting the data with a single Lorentzian line. Therefore these results are quite robust. 



We note that, for all samples, the variation of Toi(B) as a function of yB (top panels of 
Fig. [2]) displays with respect to a smooth variation, an extra bump. We have decomposed 
this variation into two parts, 701 which is quadratic in y/~B and an extra contribution 
displayed by stars in the top panels of Fig. |2J We will assign this latter contribution to the 
i^-phonon interaction contribution. We have taken for 710 a quadratic dependence of 
the form: 701 (B) = 70 + 71 x \JB + 72 x B. Though it is possible to fit F ol (B) for B lower 
than 16 T with 7 2 = 0, the resulting decomposition gives, for the i^-phonon contribution, 
a variation which does not decreases at high fields (or high energies) which is not physical. 
We are then lead to introduce this 5-dependence which, in turn, will determine the way 
the interaction decreases at high energies. This is of course dependent of the model used 
to describe the interaction as discussed in section V-B. We have used the same parameter 
72 = 0.24 for all samples which in turn gives a value for 71 = 0.55 ± 0.05 (the error bars 
covering the coefficients for all samples). This dependence on B is not explained by existing 



scattering mechanisms 



2]. Coming back to the i^-phonon interaction contribution, we find, 



for all samples a similar variation: after a threshold around 17 T, it reaches a maximum of 
2.1 ± O.lmeV (the error bars covering the values found for all samples). 

In the bottom panels of Fig. [2] (and Fig lb of the main paper), the variation of Eqi/B 1 / 2 
versus B 1 ^ 2 is displayed. These variations show a downward kink at a value of B corre- 
sponding to the onset of the i^-phonon contribution. It is important to note that, within 
the experimental errors, the i^-phonon contribution to Tq\(B) and the variation of Eqi/B 1 ^ 2 , 
plotted as a function of the energy are indeed Kramers-Kronig related. 

Though the present analysis did not focus on the higher energy transitions, we have fit 
as well over a large energy range the transitions between the n = —2 to n — 1 and n — — 1 



to n = 2 Landau levels corresponding to the energy e_i )2 = e_2,i = Vp\/2ehB x (\/2 + 1). 
For all samples this transition can be fitted with vf = 1.025 x 10 6 ms _1 . It is interesting 
to note that the width of these transitions T_ 12 (B) is found to vary linearly as a function 
of B l l 2 like r_ 12 (-B) = f]o + f]\ x \f~B without, within the experimental errors, a significant 
contribution proportional to B. Therefore it seems that the dependence on B adopted for 
7oi(-B) is specific of the n = LL. 

III. FITTING THE EXPERIMENTAL TRANSMISSION CURVES 

Up to this point the treatment of the transmission curves has been done assuming that all 
optically active layers have the same broadening parameter Tiq(B). In reality the variation 
of E m / B 1 / 2 displayed in the bottom panels of Fig. [2] lead us to conclude that this if-phonon 
interaction, reminiscent of the Frohlich interaction observed in polar semiconductors 3|, |J] , 
has to be treated as a function of the energy. In other words, for each layer simulated 
with Eq.l of the main text, the denominator hu — Eqi(B) + iTqi(B) has to be replaced 
by Tko — eoi(-B) + ry i(-B) + H(ui,B). Therefore the new pole occurs at an energy e i(B) — 



Re(E(u>,B)) with a width j i(B) + Im{T,{uj,B)). Here E(co,B)) is the self energy of the 
phonon interaction and eoi(-B) = VF\/2efiB the non re-normalized value of the transition 
energy. The model for treating the self energy will be detailed in section V: the imaginary 
part of the self energy depends on a parameter R (see Eq JT7|) . the phonon energy ftw K and 
an amplitude squared A 2 governing the electron-phonon interaction (see section V-B). For 
each sample we have used the background broadening parameter 701 (B) as obtained from the 
preceding fitting (which is sample dependent) but the parameters for the phonon-interaction 
are the same for all samples: R = 3, TujJk = 151 meV and A 2 = 1.94 meV 2 x B. The same 
Fermi velocity vf = 1.012 x 10 6 ms _1 has been taken for all samples to describe the transition 
eoi(-B). It is a little bit smaller than the value found for the transition e_i2(5) which may 



be a sign of electron-electron interactions 



!]• 



The results from the model, presented in the next sections, are compared to the exper- 
imental data in Fig. [31 Including the experimental errors, the parameter R is determined 
within 15 per cent. However the fitting is done assuming that the valley splitting Ay is 
larger than the spin splitting A5, a result which is clearly beyond the experimental uncer- 
tainties as evidenced by the comparison of both simulations with the experimental spectra 
(see Fig. 4 of the main text). As a matter of fact, our model predicts, see section IV-B, that 
the i^-phonon interaction should vanish when the filling factor v goes to zero if Ay < A$- 



Despite the different approximations made to treat the spectra, we obtain a fit of trans- 
mission curves which is acceptable with the same basic parameters (t>o, Vp, R, hu)K, A 2 ) for 
all samples. This makes us confident in the basic physical picture and the process we have 
adopted to treat these spectra. 
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FIG. 3. Comparison of experimental data (open black dots) with the simulated traces (red lines) 
using the proposed K-phonon interaction model for different values of the magnetic field. Theo- 
retical spectra (from SI-IV) were calculated with parameters described in the text, (a) sample S2 
(scale multiplied by 2), (b) sample S3 and (c) sample S4. 



IV. DERIVATION OF THE MODEL HAMILTON! AN 



A. Electron-phonon interaction 



Expressions for the electron- .fT-phonon interaction in graphene in the plane wave basis 
can be found in Ref . 5] . After performing a change of basis to the Landau level (LL) basis, 
the electron-phonon interaction can be written as: 



H e ~ ph = £ 9^c] +li<k/Cj>n , k (al ? +a^. (1) 

jnn' 
kk'g 

Here, ct n k is the electron creation operator in valley j = 0, 1 (corresponding to the K 
or K' valley), for the LL n and wavevector k in the Landau gauge, and a}- - is the K or K' 
phonon creation operator of wavevector q relative to the K or K' point. The addition of 
the valley index j + 1 is understood to be modulo 2. The matrix elements are: 
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(2) 



with 



7 <n(0 =g) 1/2 e-^ (^±i^) n '" n L«'-«(^l) whenn'>n>0 



(3) 



Jn'-i(q) = 

In these equations, 6 = a/v^ is the C-C bond length, a = 2.46 A is the lattice period, 7 = 
\^3at /2, and t = 2.7 eV is the nearest neighbor hopping parameter. (3 = —(b/^)(d / y/db) 
is a dimensionless coupling parameter. N c is the number of unit cells, M is the mass of a 



carbon atom, 00 k is the if-phonon frequency and Ib = \fh] (eB) is the magnetic length. 
The constant C takes the value of 1/2 if |n| > and \n'\ > 0, and l/y/2 otherwise. L^(u) 
denotes the Laguerre polynomial and sgn(n) = 1 when n > 1 and sgn(n) = —1 otherwise. 

B. Choice of reduced Hilbert space 

The cyclotron transitions of interest are (see Fig. 2 of main text): 

\ViU, k)) = ct i fcCj.o.fclO) 

(4) 

\mU> k )) = c J,o,fc c J,-i,fcl°>- 

Here, |0) is the ground state. We now consider candidate states |£) to be included in the 
model, such that the experimental spectra could be explained by their interactions with the 
cyclotron transition states. Such states have to satisfy two conditions: (i) they can couple 



directly to (772) via the electron-phonon interaction Eq. [T] (i.e. (£|i? e ^771,2) 7^ 0), 
and (ii) in order to describe the resonance at E = Jvjjk the states |£) must have energy 

m\o « hw K . 

Based on the form of the electron-phonon interaction, it is evident that |£) must contain 
a i^-phonon as well as intervalley electronic transitions. Condition (ii) implies that these 
transitions must be between LL of the same LL index because the electronic part of |£) must 
have the same energy as |0). We therefore consider the states 

Note that the single phonon states aj^O) satisfy condition (ii) but not condition (i). In 
Eq. El the inter-valley transitions are between the n = LL because this is the only partially 
empty LL in the experiment. 

Another class of states that can couple to the cyclotron transition states is 
{^l£<?(i> k)), &)))•■ ■}> which contain one or more acoustic phonons in addition to 

|£<f(j, k)). Here, W , create acoustic phonons. 

C. Integration over the phonon continuum 

Because of the dispersion of optical phonons at the K point [6|, and the coupling to 
acoustic phonons, the energies of the states that can couple to |r/ 2 ) form a continuum, 
bounded below by u>k- In this section, we perform an integration over this continuum. We 
start with a Hamiltonian of a general form that describes two degenerate states interacting 
with a continuum of states. For a given (j, k), 
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/ 



#1 



Eoi ui(uj k ,oj k ) ui(uj k + 5, u K ) ui(u K + 25,u K ) ■ 

E i u 2 (ujk,uk) u 2 (u k + S,ojk) u 2 (u k + 25,u k ) ■ 

ui(uj k ,uj k ) u 2 (u K ,u K ) U K V V 



Ui((jJ K + 5,U3 K ) u 2 {u K + 6,U K ) 

ui(u K + 25, oj k ) u 2 (u k + 25,u k ) 



\ 



uj k + 5 



ui K + 28 



(6) 



The first two columns (rows) represent the cyclotron transition states |?7i), \rj 2 ). The 
rest of the columns (rows) represent the continuum of states that contain phonons. Here, 
5 > is a small positive energy. The function Ui(u,UK) (may be obtained from Eqs. [fl [2]) 
describes the interaction between \r]i) and a |£) state of the continuum of energy u, while v 
describes the interaction within the continuum. 



We now treat the problem in two successive approximations. The first one assumes that 
the phonons have no dispersion, which allows us to introduce an effective coupling constant V 
and work in a restricted Hilbert space having effectively one phonon. The second assumption 
takes into account the dispersion of phonons as a correction to the first approximation. This 
is done by introducing a self energy term in the Green's function (see next section for the 
treatment of the Green's function). 



An effective Hamiltonian in our restricted Hilbert space for dispersion-less phonons can 
be written as: 
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H 



^01 



£ i V 2 



(7) 



\v; v 2 * u K J 

where the values of V\ and V2 are chosen such that H and H\ have the same eigenvalues 



when 5 = 0. We have found that V\ = — V 2 = V = w (see below for explicit 
expressions of this parameter.) 

We then introduce the effects of the phonon continuum using second-order Lowdin pertur- 
bation theory |7|. This is done by introducing a self energy H ph (uj) in our Green's function 
treatment of the optical response functions in the restricted Hilbert space of if 8|. The 
condition imposed by the model is that the Green's functions (ou — H — Y, p h(uj) — i0 + )~ 1 and 
(cj — Hi — i0 + )~ 1 are equal in the space spanned by j^)- Explicit expressions for the 
different components of Ej (u) are given in the next section. 

In the definition of the parameter V, the prime (') on the sum indicates that not all q 
vectors have to be summed over, because only a fraction < / < 1 of the total number of 
states are available, depending on the filling factor. Then, 



all q 



(8) 



The constant A 2 = Y^aii qdiqdiq = 4^3mL^ can ^ e obtained from Eq.JEJ In reality this 
has to be corrected by the introduction of electron-electron interactions [6j. To obtain the 
value of /, we examine the occupation numbers of the n = LL of the K and K' valleys. 

Let us define v to be the filling factor, including spin and orbital degeneracies. For 
example, v = at charge neutrality and v = 2 if the carrier concentration is such that the 
LL n = is fully occupied. For the case of Ay > A s , filling factor v < 1 (Fig. H^), we 
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FIG. 4. a) Schematic diagram of the Landau levels when Ay > Ag and v < 1; b) Schematic 
diagram of the Landau levels when Ay > Ag and 1 < v < 2 and c) Schematic diagram of the 
Landau levels when Ay < Ag and v < 1, 



have v7 = 1 f° r t ne spin up transitions, and y/J = y/1 — v for the spin down transitions. 



Whereas for Ay > As, filling factor 1 < v < 2 (Fig. Hb), we have vT = a/2 — v for the spin 
up transitions, and \fj = for the spin down transitions. 

On the other hand, if Ay < As, the expressions are the same as the case of Ay < As 
when the filling factor is 1 < v < 2, but when v < 1: we have y/J = \fv for the spin up 
transitions, and y/J = for the spin down transitions (Fig. Hfc). Therefore, in that case, for 
undoped graphene layer, the interaction disappears. 



V. EVALUATION OF OPTICAL CONDUCTIVITY 
A. Green function formalism 

The Green function formalism for the optical conductivity Reajj ~ Im (0\M^GMj\0) has 
been introduced in the main text. Here, we provide some explicit expressions for the form 
of the optical conductivity. In order to facilitate the discussion of the self energy S p / l (o;), we 
express the Green's function matrix in the basis of eigenstates of the model Hamiltonian H 
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(Eq. ED. 



^G' n ^ 



G' 



G' 22 



V 



, Ga 



1 



J J 



LO — Sj — IS 



8 = 0' 



(9) 



G' 33 J 

The primed (') quantities are expressed in the basis of the eigenstates of H. The eigen- 
values Ej of the Hamiltonian H (Eq. UJ) are 



£i = -Em 

£2,3 = \ (E 01 + ^ ± V /(E 01 - huj K ) 2 + 8V 2 ) 



(10) 



The Green's function matrix (G) in the original basis (i.e., the same basis as Eq. [7J can be 
recovered by a unitary transformation G = UG'U^ 1 . 



U 



(l 


-V 




-V 




V2 


A /2V 2 +( J B i- 


-e 2 ) 2 


y/2V 2 +(E i- 


s 3 ) 2 


1 


V 




V 






y/2V 2 +(Eoi- 


-e 2 ) 2 


y/2V 2 +(Eoi- 


-e 3 ) 2 




Eoi-62 




Boi-£3 




^2V 2 +(E 01 - 


S2? 


y/2V 2 +(Eoi- 


-es) 2 



(11) 



Then, for instance, the a xx component of the optical conductivity is related to (0\M}.GM x \0) 
(0|M,J[/G / [/~ 1 M r |0), where M\ = (M l2 ., M 2x , 0) is a row vector containing the optical matrix 
elements. 



B. Discussion of S 



ph 



We introduce self energy terms into the diagonal elements of G' 



jj 



oo — ej — i0 + lo — Ej — £j(o;) 
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(12) 



From second-order Lowdin perturbation theory (see previous section), we have found 



ImSi = 



Im E 2 = ojk) 



Im E 3 = wx) 



where the imaginary part of the self energy has the form 



ImE(w) = 9(co — ojk)(u(oj, lojc))" 



(13) 



(14) 



Here, u(uj,ujk) = \ui(w,wk)\ = ^(w, cok)\- In the basis of Eq. [7J the Green's function can 
be written as 



G 



Gn — G 



1 



22 



t 1 f 1 ^ 
Uii Ori2 Cri3 



(jr21 G*22 G23 
^31 G32 G 3 3 / 
1 



1 



G 



2\(jj — E 01 -i0+ u- E m - i0+ - E(w) 
I / 1 1 
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(15) 



(16) 



2\ lu- E 01 -iO+ lu- E QX ~ i0 + ~ 
Im S(w) has a threshold at uk that arises irrespective of the form of the coupling matrix 

element u. In this model, the threshold is always located at the lower bound of the energies 

of the continuum of states. We have chosen, motivated by the experimental data (see top 

panels of Fig. [2]), for u 2 the form 



(u(lo, cor)) 2 = Dexp(—R(u/uK — 1)) 



(17) 



Using the sum rule for Im(G) together with the constraint that Hi must reduce to H 
in the limit where R — > 00 (the limit where the continuum collapses to a single energy 
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FIG. 5. Top panel: Imaginary part of the self-energy (£2 or £3 in Eq. IT3|) with v = 0. Bottom 
panel: Real part of the self-energy with v = as deduced from Kramers-Kronig transformation 
of the Imaginary part. The behavior of the real part of the self-energy, as a function of energy, is 
quite similar to that of Eqi/B 1 / 2 when adding the contribution of the non re-normalized value of 
the transition energy €oi(B)/y/~B = vpV^eh = 36.8meVyT 1//2 . (See bottom panels of Fig. [2] and 
Fig. lb of the main text.) 



ujk)i the value of D can be fixed at D = (ttRV 2 ) / '(hux). The model is not sensitive 
to the exact functional form of Eq. [0] as long as the qualitative features such as the 
exponential decay and the presence of a threshold at uj — u>k are preserved. We then have 
two independent parameters to describe the interaction, V 2 and R. As derived in Eq. [HI 
V 2 = A 2 f and we will use the value A 2 ( meV 2 ) = 1.94 x B(T) [6J. The energy range 
of u(u,Uk) contains contribution from the i^-phonon dispersion and multi-acoustic phonon 
processes. The only remaining fitting parameter is R. In the interpretation of the continuum 
as states containing multiple acoustic phonons, the exponential factor in Eq. [17] describes the 
decreasing probability of multiple acoustic phonon emission far from the threshold energy. 
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With this simplified model, and knowing the experimental results, we have chosen the value 
R = 3 as sketched in Fig. [5] which in turn is coherent with the value of the linear B term 
used for 701 (B). This also gives reasonable values for ReS(cj) (bottom panel Fig. E]) as 
compared to the experimental values. Of course the extension in energy corresponding to 
this R value is quite high but it does not play a role in the actual fitting process limited 
to about 200 meV (see Fig. [3]). Though, within the experimental errors, the fit of data is 
acceptable with R — 3, it could be also accepted for R values, as well as the 72 values of 
7oi(-B) differing by 15 per cent. Knowing that the spin splitting A$/B ~ 0.12 — 0.15 meV/T 
and that the the introduction of Ay oc B should result in an extra broadening of the E i 
transition, the introduction of the 72-term is explained naturally by the fact that it should 
be larger than A$/B. 

We note finally that the introduction of the interference term G 12 (Eg JIp 7 ]) (which corre- 
sponds to interference effects between the cyclotron transitions at the two valleys) is essential 
to reproduce the experimental results. 
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